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As shown by Taubes, in the Bogomol'nyi-Prasad-Sommerfield limit the SU(2) Yang-Mills- 
Higgs model possesses smooth finite energy solutions, which do not satisfy the first order Bo- 
gomol'nyi equations. We construct numerically such a non-Bogomol'nyi solution, corresponding 
to a monopole-antimonopole pair, and extend the construction to finite Higgs potential. 
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I. INTRODUCTION 



SU(2) Yang-Mills-Higgs theory, with the Higgs field in the adjoint representation, possesses magnetic monopole 
and multimonopole solutions. The solutions with unit magnetic charge are spherically symmetric |l] [|]. In 
contrast, multimonopole solutions cannot be spherically symmetric |5| and possess at most axial symmetry |^-^]. 
In particular, for magnetic charge greater than two, solutions with no rotational symmetry exist |Tc| ]. 

In the limit of vanishing Higgs potential, the BogomoPnyi-Prasad-Sommcrficld (BPS) limit, monopole and mul- 
timonopole solutions satisfy the first order Bogomol'nyi equations [[ll} as well as the second order field equations. 
They have minimal energies, saturating precisely the Bogomol'nyi bound. 

In the BPS limit, monopole and axially symmetric multimonopole solutions are known exactly In 
contrast, for finite Higgs potential monopole and axially symmetric multimonopole [fl2| solutions are known 
only numerically. But even in the BPS limit, numerical construction of axial multimonopole solutions Q pre- 
ceeded their exact construction , and multimonopole solutions without rotational symmetry are only known 
numerically [fiof . 

As shown by Taubes |L^|, "there is a smooth, finite action solution to the SU(2) Yang-Mills Higgs equations 
in the Bogomornyi-Prasad-Sommcrficld limit, which does not satisfy the first-order Bogomol'nyi equations". We 
here construct numerically such a non-Bogomol'nyi BPS solution, first found in |Q. This solution possesses 
axial symmetry and corresponds to a monopole-antimonopole pair. We extend the construction to finite Higgs 
potential. 

We review the SU(2) Yang-Mills-Higgs model in section II and the axially symmetric ansatz for the monopole- 
antimonopole solution in section III. We analyze the magnetic charge of the solution in section IV and present 
numerical results in section V. In section VI we present the conclusions. 

II. St/ (2) YANG-MILLS-HIGGS MODEL 

We consider the SU (2) Yang-Mills-Higgs Lagrangian 

-£ = j ' ^Tr(F^F^) + \Tr(D^D^) + ^Tr((<£ 2 -r? 2 ) 2 )}^ (1) 

with field strength tensor of the su(2) gauge potential — ^T a A^, 

F tllJ =d ll A v -d„A tJ , + i[A tJ „A v ] , (2) 

and covariant derivative of the Higgs field $ = T a <fi a in the adjoint representation 

D^ = d^ + i[A„,^} , (3) 

and g denotes the gauge coupling constant, A the strength of the Higgs potential and rj the vacuum expectation 
value of the Higgs field. 

The Lagrangian (Q) is invariant under local SU{2) gauge transformations g, 

Af, — > g^g -1 + idf+gg- 1 , 
$ — > g$g _1 . 

(4) 

Static finite energy configurations can be characterized by an integer topological charge Q 

Q= 4^J Tr {F,3i}t$} ' (5) 

corresponding to the magnetic charge m = Q/g. In the BPS limit the energy E of configurations with topological 
charge Q is bounded from below 
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E>^9_. (6) 

Monopole and multimonopole solutions satisfying the Bogomol'nyi condition 

F« = e tjk D k r (7) 

precisely saturate the lower bound (||). 

Here we construct a solution, which corresponds to a monopole- antimonopolc configuration and therefore carries 
Q = 0. It has finite energy E > 0, and thus, in the BPS limit, corresponds to a non-Bogomol'nyi solution of the 
SU(2) Yang-Mills-Higgs field equations. 

III. STATIC AXIALLY SYMMETRIC Q = ANSATZ 

We choose the static, axial ly sym metric, purely magnetic Ansatz employed in for the monopole- 
antimonopole solution and in |l5|,[16| for the sphaleron-antisphaleron solution of the Weinberg-Salam model. 
Here the gauge field is parametrized by 

Ao = 0, A r = ^r V) A e = { l^Mr 9) A v = (h^ + (1 - HA rf A , (8) 

2gr g g \ "J 

and the Higgs field by 

$ = rj(*iT< 2 ) +$ 2 r e (2) ) . (9) 

(2) (2) 

The su(2) matrices 7V , Tq and t v are defined in terms of the Pauli matrices ti,T2,T3 by 



= sin 26>(cos ipTi + simpT2) + cos29t 



( 2^ 

Tq ' = cos 2#(cos ipTi + sin (pT 2 ) — sin20r 3 , 
t v = -sin^ri + cos (fiT 2 , (10) 

and for later convenience we define 

T p = COS ifTl + Sin LfT 2 ■ (11) 

We change to dimensionless coordinates and Higgs field by rescaling r — » r/(grj) and <& — > 77$, respectively. 
Then this Ansatz leads to the field strength tensor 

F re = -^-(dgH 1 + 2rd r H 2 )T v , 
2r 

F rv = — ((sin 20ffi - 2sin0.ff x (l - HA -2sm0rd r HA r< 2) 
2r I 

+ (cos 20i? 1 + 2 sin 0i? i/f 3 + 2 sin 6rd r H A ) r e (2) | , 
F flv = -- {(2 sin 26(H 2 - 1) + 2 cos 6>iJ 3 - 2sin6»iJ 2 (l - iJ 4 ) + 2 sm6d e H 3 ) 

+ (2cos2e(H 2 -l)+2cos6{l- HA. +2smdH 2 H 3 - 2sin0d e HA.T e :2) \ , (12) 
and the covariant derivative of the Higgs field 

D r $ = - {(r3 r $i + fli* 2 ) ^ (2) + (rd r $ 2 - ffi$i) r e (2) } , 
D $ = (a e $i - 2ff 2 $ 2 ) r r (2) + (9 e $ 2 + 2tf 2 $i) r e (2) , 

D v $ = {(sin2(9-2sin6>(l-F 4 ))$i + (cos 29 + 2 sin 6H 3 ) $2} r v . (13) 
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The dimcnsionless energy density then becomes 



- I r <! — F 2 g + - . 2 F 2 + - . 2 Fg 
r z r 2 sin 6 Y r 4 sin 6* * 



-Tr |(A.$) 2 + 4j(A^) 2 + — ^(^*) 2 } + A((|$ I 2 - I) 2 , 



(14) 



where |$| = y / $f + $f denotes the modulus of the Higgs field. 

For finite energy configurations the modulus of the Higgs field has to be one at infinity, whereas the covariant 
derivatives of the Higgs field have to vanish at infinity. These conditions lead to jl4|] 

r — > oo : $i — ► 1 , $ 2 > , (15) 



Hi 



0, H 2 ^0 7 H 3 



suit 



1 - Ha 



cost 



Substituting the asymptotic expressions for the gauge field functions into the Ansatz 
potential approaches a pure gauge at infinity 



(16) 

shows, that the gauge 



0, 

— sin6*(cos 0T n 



= -i<9 v gg f 



(17) 



where g = expji^T^} rotates the Higgs field at infinity to a constant, (g'&g^)^ = T3 p"7| . Reexpressing the 
topological charge (g) as a surface integral, we find Q = for configurations obeying (|15|), ([l6|), provided, the 
configurations are sufficiently regular. 

Inserting the Ansatz (^|), (^J) into the general variational equations leads to a system of six coupled non- linear 
partial differential equations for the four gauge field functions H t and the two Higgs field functions <fv The 
same system of partial differential equations is obtained by inserting the Ansatz directly into the Lagrangian and 
calculating the variation with respect to the functions Hi and showing that the Ansatz (||), (||) is self-consistent. 

The Ansatz is not a priori well defined on the z-axis and at the origin. However, for solutions of the field 
equations, we have performed an expansion of the gauge and Higgs field functions near the z-axis and near the 
origin |L8|]. Inserting these expansions into the Ansatz we find that the gauge potential and the Higgs field are 
well defined and (at least) twice differentiable on the z-axis and at the origin. 

At the origin we find 



Ay 

A z 



93 

-~2 y Tl 



|(. 9l + 2 52 ) + |(2 :c 2 



Z -(g 1 + 2g 2 ) + 9 -?-(2y 2 
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(.92 +g3z)(yn -xt 2 
(3p 2 - 7: 



z 2 ) 



z 2 ) 

93 



T2 ' 

xyr 2 



9iyzTz , 
-54XZT3 , 



10 



96 (p 



4z' 

00 35 



(xti + yr 2 ) 

-g 6 )(3p 2 -2z 2 )+g 7 (p 2 -2z 2 ) 



T3 



where <fio,gi are constants. Therefore the Ansatz allows for a non- vanishing Higgs field, <i>(r = 0) 
origin. 

Near the z-axis the gauge field functions behave like 



(18) 

6 r 3 , at the 



Hi = hu(r) sin (9 - 
H3 = h 31 (r) sin9 



H 2 = f{r) + h 22 (r) sin 2 6 + .. 
H A = f{r) + h A2 {r) sin 2 6 + .. 



(19) 



while the Higgs field functions behave like 
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$1 = 4>{r) + cj)i 2 {r) sin 2 6 + . . . , $ 2 = 021 (r) sin6» + . . . , 



(20) 



where . . . indicate higher order terms in sin6>. At the nodes zq of </>(r), the modulus of the Higgs field vanishes. 
Therefore, these nodes correspond to the locations of monopoles. 
The Euler-Lagrange equations possess the discrete symmetry 

z ^ -z , ft -> -H lt H 2 -> H 2 , H 3 ^ H 3 , 1 - H 4 -> -(1 - H 4 ), $ x -> * ls $ 2 -> -$ 2 . (21) 

When the solutions possess the same symmetry, we expect for each node zo on the positive z-axis a node — zo 
on the negative z-axis, i. e. the monopoles come in pairs. For a solution with total magnetic charge Q — 0, then 
half of the monopoles carry negative magnetic charge and half of them positive magnetic charge. In the simplest 
non-trivial case the solution then describes a monopole-antimonopole pair. 

For such a monopole-antimonopole pair we expect a magnetic dipole field for the asymptotic gauge potential. 
Deriving the asymptotic behaviour of the gauge potential, we find that the gauge field function H 3 decays like 
0(1/ r) at infinity, while the other gauge field functions decay exponentionally. In particular, in the gauge where 



the Higgs field is asymptotically constant |17| i. e. $ — > r 3 , 



ff 3 = — sin0, (22) 
r 



which leads to the asymptotic gauge potential 



Ai = C m &^fiin, (23) 

representing indeed a magnetic dipole field. 

IV. MAGNETIC CHARGES OF THE Q = CONFIGURATION 

As a consequence of the vanishing topological number, Q = 0, the configuration we are considering carries 
zero net magnetic charge, m = 0. In the following we demonstrate, that this field configuration still possesses 
magnetic charges. 

Let us parameterize the Higgs field as 

$ = $it p + $ 2 r 3 , (24) 

with 

$1 = sin20$i +cos20$ 2 , $2 = cos2#$i - sin26»$ 2 , (25) 
and define the normalized Higgs field by 

$ = cosQTp + sin ar 3 , (26) 

with 



$i = |$| cosa , $ 2 = |<I>| sina , |$| = \J $f + $| = J$l + $| . (27) 

Thus, <3E> maps any closed surface S in coordinate space to a 2D sphere in isospin space. We define the degree of 
the map as 



2V(S) 

where V(S) is the volume of the surface S. 



piS) = wfs)L Tr {* d * Ad *}> (28) 

Let us first calculate the degree of the map of a 2D sphere S 2 of radius r centered at the origin. We find 
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P iS2) (r) 



1 

4vT 



27T f 7T 



dg sin adOdp 



• sin a 



(29) 



We now anticipate that the modulus of the Higgs field possesses two zeros located at zq and — zq on the positive 
and negative z-axis, respectively. Then the function a possesses discontinuities on the z-axis at zq and — Zq, 
i. e. a = —3tt/2 for — oo < z < —zq, a = —n/2 for — zq < z < zq and a = tt/2 for zq < z < 00. For any r we find 
then p( s ")(r) = 0, thus the map S 2 — ► S 2 can be contracted to the trivial map. 
Next we consider the integral over the half sphere 



Hli={(r,0,<p) I r fixed ,0 e [0, |] )¥ >e [0,2tt]} 



Note that the boundary of is a circle in the icy-plane, 



I 



r fixed , = — , ip s 



[0,2tt]} 



(30) 



(31) 



On the boundary the map $ is constant, $(r, tt/2, ip) — — t 3 . We now compactify ff? to a 2D sphere 5?, by 
identifying the boundary of with the south pole of . Calculating the degree of the map we find 



(r) = -TZ 



■in 



2tt r ir/ 2 



dg sin adOdtp 



— sin a 
2 



0=7r/2 



)=0 



0(r - «b) • 



(32) 



Thus, the degree of the map vanishes if the location of the zero z Q is not inside the sphere Sl_ and equals one 
otherwise. An analogous calculation for the lower half sphere H 2 , leads to 



\r)=-e(r-z ) 



(33) 



Let us now compare with the magnetic charges m in the upper and lower half spaces. To this end, we consider 
the 't Hooft electromagnetic field strength tensor 



Tu. v = Tr QFilv 



or equivalently 



The magnetic charge inside a closed surface S can be expressed as 

*F u dx dx 



V(S) 



(34) 



(35) 



(36) 



Note, that the integration of the first term in Eq. (pSf) leads to the degree of the map. 

For the upper half space we define the closed surface S+ = H\ U D 2 , where D 2 is the disk of radius r in the 
xy-plane centered at the origin. Taking into account the correct orientation of the surface element, we find 



,(5+) 



1 




V(S+) 




= + 


-1 

4tt 


= p^ + 


— 

4tt 



Te^dQdp 



J- r <pdrd<p 



D 2 



Tr 



H 



+ 
2tt 



{dofiAyijdBdkp - J 2 Tr {d r ,(i>A v )}dr'd<p 



Tr 



{(i>A v ) \ r } 



9=0 



dip 



Tr{($A v ) | 9=f } 



dtp 



(37) 
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where we have used that Tr{$Ag} — Tr{QA r } = 0, and that the first term in Eq. ( pS5| ) does not contribute to 
the integration over the disk D 2 . From symmetry considerations we know that A^lg^o = 0, A ip \g = * = H3T3 and 
# 3 (r = 0) = 0. Hence we find (see also Appendix A) 

m (s+) = p {H%) + J_ ^27r(2sina J ff 3 )| e=f - 27r(2sinaff 3 )| e= , 

= pW> . (38) 
Analogously, for the lower half space we find 

These calculations show indeed, that the configuration possesses two magnetic charges with opposite sign, 
located on the positive and negative z-axis, respectively. In general, for closed surfaces which contain only one of 
the locations of the zeros of the Higgs field, integration of the field strength tensor normal to the surface yields a 
non- vanishing magnetic charge. In constrast, for a surface enclosing both charges, their contributions compensate, 
yielding zero net magnetic charge. 

V. NUMERICAL RESULTS 

The Ansatz Eqs. @ is form invariant under abelian gauge transformations g = e %T l 2T v , with the gauge 
and Higgs field functions transforming as 

-^^-^ + d r r, 

r r 
2# 2 -> 2# 2 + d g T , 

H 3 ~^) ^cosr(# 3 -^) + sinr(l-# 4 -cos0) , 
2 sin 9 J \ 2 sin 9 ) 

( cos 20 

(1 - #4 - cos 0) -> cos r (1 - #4 - cos 0) - sinT #3 

\ 2 sin / 

$1 — > cos T$i + sin T<i>2 , 

$2 — > cosr$ 2 - sinT^! . (40) 
To find a unique solution we have to fix the gauge and choose the condition 

G t =-j (rdrHi - 2dgH 2 ) = . (41) 

The set of partial differential equations is then obtained from the Lagrangian Eq. (|l|) with the gauge fixing term 
£Gf added, where £ is a Lagrange multiplier. 

This set of partial differential equations is solved numerically subject to the following boundary conditions, 
which respect finite energy and finite energy density conditions as well as regularity and symmetry requirements. 
These boundary conditions are at the origin 

#i(O,0) = #3(0,0) = 0, #2(0,0) = # 4 (O,0) = 1 , (42) 

sin20$i(O,0) + cos20$ 2 (O,0) = , <9 r (cos20$i(O, 0) - sin20$ 2 (O, 0)) = , (43) 

at infinity 

#i(oo,0) = #2(00,0) = , #3(00,0) = sin , (1 - #4(00,0)) = cos (44) 

$1(00, 0) = 1, $ 2 (oo,0)=O, (45) 
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and on the z-axis 

Hi(r,6 = 0,ir) =H 3 (r,6 = 0,tv) = d e H 2 (r,8 = 0,n) = d e H 4 (r,0 = 0,n) = , (46) 



$ 2 (r, = 0,tt) = 9 e $i(r, (9 = 0, tt) = . (47) 

We have constructed monopole-antimonopole solutions for a large range of values of the Higgs coupling con- 
stant A. The numerical calculations were performed with the software package CADSOL/FIDISOL, based on 
the Newton- Raphson method fl9| . For vanishing Higgs coupling constant the monopole-antimonopole solution 
corresponds to a non-Bogomol'nyi BPS solution, for which our results are in good agreement with those of 
Ref. @. 

In Table 1 we present the normalized energy of the solutions E/Anr) for several values of A and compare with 
the energy Ei n f /Airrj of a monopole-antimonopol pair with infinite separation, corresponding to twice the energy 
of a charge-1 monopole. For all values of A in Table 1 the energy of the monopole-antimonopole solution is less 
than the energy of a monopole-antimonopole pair with infinite separation. 



A 


E/4nr] 


E iDi /4wri 


d 


4>o 


Cm 





1.697 


2.000 


4.23 


0.328 


2.36 


0.001 


1.830 


2.053 


3.48 


0.381 


2.07 


0.01 


2.015 


2.204 


3.34 


0.489 


1.84 


0.1 


2.330 


2.498 


3.26 


0.791 


1.71 


0.2 


2.442 


2.613 


3.24 


0.886 


1.69 


0.5 


2.596 


2.776 


3.11 


0.961 


1.62 


1.0 


2.713 


2.900 


3.0 


0.986 


1.57 


10.0 


3.042 


3.241 


3.0 


0.9996 


1.55 



Table 1 

The energy of the monopole-antimonopole solution as well as the energy of two infinitely separated Q = ±1 
monopoles, the distance d between the locations of the monopole and antimonopole, the modulus of the Higgs 
field at the origin, <^o, and the dimensionless dipole moment C m are given for several values of the Higgs coupling 
constant A. 



In Fig. 1 we exhibit the modulus of the Higgs field \&(p, z)\ as a function of the coordinates p — x 2 + y 2 and 
z for A = and A = 1. The zeros of |$(p, z)| are located on the positive and negative z-axis at ±zo ks 2.1 for 
A = and at ±zo ~ 1-5 for A = 1. The distance d of the two zeros of the Higgs field decreases monotonically 
with increasing A, as seen in Table 1. 

Asymptotically z)| approaches the value one. For A > the decay of the Higgs field is exponentially. 
The value of the modulus of the Higgs field at the origin increases monotonically with increasing A (see Table 1). 
While <j>o — 0.328 for A = 0, <j>o is already close to one for A = 1. In the limit A — > oo we expect the modulus 
of the Higgs field to be equal to one everywhere, except for two singular points on the z-axis, representing the 
locations of the monopole and antimonopole. In contrast, the angle a should remain a nontrivial function in this 
limit. This would then be similar to the result found in for the charge-2 multimonopole. 
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Fig. 1 (a) Fig. 1 (b) 




z 



FIG. 1. The modulus of the Higgs field as a function of p and z for A = (a) and A = 1 (b) 

In Fig. 2 we show the energy density of the monopole-antimonopolc solution as a function of the coordinates 
p = \J x 2 + y 2 and z for A = and A = 1. At the locations of the zeros of the Higgs field the energy density 
possesses maxima. For A = 1 the maxima are more pronounced compared to the case of vanishing A. At large 
distances from the origin the energy density vanishes like 0(r~ 6 ). For intermediate distances from the origin the 
shape of equal energy density surfaces looks like a dumb-bell. For smaller distances the dumb-bell splits into two 
surfaces. 

Near the locations of the zeros of the Higgs field the equal energy density surfaces assume a shape close to a 
sphere, centered at the location of the respective zero. This presents further support for the conclusion, that at 
the two zeros of the Higgs field a monopole and an antimonopole are located, which can be clearly distinguished 
from each other, and which together form a bound state. This is in contrast to the axially symmetric charge-2 
multimonopole solution, where the individual monopoles cannot be distinguished, and where, in fact, the Higgs 
field has only one (double) zero. 




z 



FIG. 2. The the dimensionless energy density as a function of p and z for A = (a) and A = 1 (b) 

Considering finally the electromagnetic properties of the monopolc-antimonopole solution, we observe, that the 
dimensionless dipole moment C m decreases monotonically with increasing A (see Table 1). 
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VI. CONCLUSIONS 



We have considered static axially symmetric solutions of the SU (2) Yang-Mills-Higgs model, residing in the 
vacuum sector. These solutions represent monopole-antimonopole pairs. The modulus of the Higgs field possesses 
zeros at the locations of the monopole and antimonopole. Clearly distinguished from each other, the monopole 
and antimonopole together form a bound state, which carries a magnetic dipole moment and net zero magnetic 
charge. However, this bound state is unstable, corresponding to a saddle point p^Jl^ |. 

We have constructed the monopole-antimonopole solutions numerically for various values of the coupling con- 
stant A, representing the strength of the Higgs potential. With increasing A, the energy E of the pair and the 
ratio E/E{ n { increase and the distance d between monopole and antimonopole as well as the magnetic dipole 
moment C m decrease, while the energy density becomes more localized around the monopole and antimonopole 
locations. 

In the BPS limit, the SU(2) monopole-antimonopole solution does not satisfy the first order Bogomol'nyi 
equations fl4|| . Hopefully, the numerical solution will be of help in constructing this non-Bogomol'nyi solution 
analytically. Recently also non-Bogomol'nyi SU(N) BPS solutions, corresponding to monopole-antimonopole 
configurations, have been found |20|. These solutions, however, are spherically symmetric. 
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VII. APPENDIX A 

At first glance it seems to be obvious, that the second term in Eq. ([35]) does not contribute to the surface 
integral, because it is the curl of the "gauge field" = Tr{<tA^}. However, <3E> is not continuous on the z-axis 
and introduces a singularity in the "gauge field" A^; its curl may contain ^-functions if the singularity is strong 
enough. 

To examine the singularity we expand the functions near the singular point f s = (0, 0, Zq), 

$1 = fo + h sin 2 , $2 = fi sin 9 , f = c a (z - z ) 
$! = sin0(/i + 2/ ) , $ 2 = f + sin 2 0(/ 2 - 2(/ + /i)) 
H 3 = sin 9H 31 , 1 - Hi = H iQ + sin 2 0H A2 , 



where Co, /i, f2, H31, H40 and H42 are constants. This leads to 

. fx + 2/q . /o + sin 2 0(/ 2 -2(/o + / 1 )) 

cos a = smt/ — , , sma = 

2 1 c,;^ 2 at i / 12 1 _• 2 n -f2 



fS + saSOfi ^/ 2 + sin 2 0/ 2 

, _ 2y ^31/0 + ^40/1 2 - 2x ^3l/0 + ^40/l 



(48) 



(49) 



Z y/z*f§ + (x*+y*)(fg + f?) " z y/z*fg + (x*+y*)(f* + f?) 

Next we define z — z — zq and keep only terms linear (quadratic under the square root) in x, y, z, 

A x ^2^^= Hiofl A y = -2^-== Hioh (50) 

zo ^/44^TJJWTW) z y/<%z*z2 + f*(x*+y*) 

Now we introduce spherical coordinates x — fsin#cosy>, y — fsin^sin^, z = fcosO centered at the singular 
point. With respect to these coordinates the components of the "gauge field" become 
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A f = 0, A 3 = 0, A^ = -2sm 2 e— Hmh (51) 

Then we find for the curl 

Fe v = --de (sin 2 S . Hioh ) . (52) 

Consequently the surface integral Fe v d6d<p over a sphere centered at the singular point vanishes. Thus 

the singularity of the "gauge field" is too weak to introduce 5-functions in the electromagnetic field strength 
tensor T^ LV . 
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